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ABSTRACT 

In most cosmological theories, the galaxy mass function at small masses is related to the 
matter power spectrum on small scales. The circular velocity function (a quantity closely related 
to the mass function) is well-studied for dwarf satellites in the Local Group. However, theoretical 
predictions and observational measurements are difficult for satellite galaxies, because of ram 
pressure and tidal stripping. By contrast, isolated dwarf galaxies are less affected by these 
processes, and almost always have enough 21cm emission to trace their dynamics robustly. Here, 
we use isolated low mass dwarf galaxies from the Sloan Digital Sky Survey (SDSS), with measured 
21cm widths, to test cold dark matter cosmology. We find consistency between the predicted and 
observed number density of isolated galaxies down to Knax 50 km s~^. Our technique yields 
a direct test of small-scale cosmology independent of the Lyman-a forest power spectrum, but 
our sample is currently statistically less powerful: warm dark matter particles heavier than 0.5 
keV cannot be ruled out. Our major systematic uncertainty is the surface brightness limit of the 
SDSS. Bhnd HI surveys, such as the ALFALFA survey on Arecibo, are expected to uncover a 
much larger number of isolated low mass galaxies, will increase the power of our constraints at 
small scales, and will propel the study of isolated galaxies to low masses previously attainable 
only in the Local Group. We use our sample to explore dwarf galaxy formation as well, finding 
that the Tully-Fisher relation for dwarf galaxies is a strong function of environment, and that the 
baryon (stellar plus neutral gas mass) fraction is only a weak function of galaxy mass. Together 
with the strong dependence of gas fraction on environment, these results indicate that for dwarf 
galaxies, gas loss and the end of star- formation are dominated by external, not internal, processes. 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — cosmology: observa- 
tions 



Introduction: the mass function and cosmology 



A critical test for any theory of cosmology and structure formation is whether it correctly predicts 
the galaxy mass fu nction. The current Cold Dark Matter model with a cosmological constant (ACDM; 



Sheth et al. 


2001 


Jenkins et al. 


2001 


Reed et al. 


2003 


— „, 

Yahaei et al. 



20041 ). A similar mass function describes the distribution of surviving subclumps withir i the dark matter 



halos ("subhalos") that most investigators associate with the sites of galaxy formation I Colfn et al. 19991 : 
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: course, 

extending this model to predict the number density of galaxies in the Universe requires including physical 
effects such as gas cooling, star-formation, supernova feedback, and possibly the formation of supermassive 
black holes and their feedback, that are too complex to follow exactly in numerical predictions. Nevertheless, 
in the last three decades numerous approximate approaches to the problem have lent us some understanding 
of what the correct predictions might be and what physical processes a successful model may involve (e.g., 



Rees & Ostriker 1977: White & Rees..l978; White & Frenk.1991: Blanton et al. 


1999i: iSomerville et al. 


2001; 


Benson et al. 2003HRobertson et allboOS: 


Croton et al.l 


2006). 



Recently, various investigators have tested the halo mass function at high masses and found consistency 
with pred ictions based on Wilkinson Microwave Anisotropy P robe three-year (WMAP3) cosmological pa- 



rameters ( Bahcall et al. 



2003 



Rines et al 



2006 



Eke et al.ll2006f ). In addition, for individual, high luminosity 



galaxies one can perform a similar test. If galaxy luminosity is related monotonically (with some scatter) to 



halo m ass, then ACDM predicts the weak-lensing signal as a function of galaxy abundanc e. iTasitsiomi et al 
(2004) have verified this prediction for galaxies with Lr > Lr^*, and Seliak et al. ( 20051 ) have performed a 
similar test extending down to about O.IL,.^*. 

Below O.lLj. the luminosity function is closer to a power law, with a slope that varies somewhat over 
luminosity but is never steeper than ab out oc L~^'^ at most, significantly shallow er than the prediction for 
the number as a function of halo mass ( Trentham et al. 2005 : Blanton et al. 2005h . If galaxies and subhalos 
are associated one-to-one, then the mass-to- light ratios of galaxies must increase substantially with decreasing 
luminosity. However, data on galaxies at the lowest luminosities has in the past been rather scarce. At circular 
velocities below ~ O.lLr, nobody has demonstrated consistency between the ACDM prediction and the 
galaxy mass func t ion. S ome valiant attempts exist, but rely on ex trapolating t he Tullv &: Fisher ( 1977) a nd 
Faber fc Jackson ( 19761 ) relations into the low luminosity regime ( Desai et al. 2004 : Goldberg et al. 2005). 



The most extreme exam ple of this issue is the "substructure problem" in the Local Group (jKlypin et al 
19991 : iKravtsov et al.ll2004b[ ). ACDM predicts hundreds of low mass satellites of the Milky Way, but only 



approximately twenty are k n own (thoug h the numb er is growing month by month; 



Ibataetal 



1995t iMatej 



1998t IWillman et"aDl2005allbl : iBelokurov et allbooel : IZucker et aLll2006alb[ ) . Isimon fc Gehal |2007l ) show that 
although the newly discovered galaxies have eased this discrepancy, a factor of 2-4 difference remains even 
when reasonable baryonic physics has been taken into account. This test of cosmology probes the lowest 
mass galaxies possible, and so is extremely sensitive to differences in the mass function slope. However, both 
predictions and observations are quite difficult in the vicinity of a luminous galaxy. From the point of view 
of theory, a number of processes, such as tidal stripping, ram pressure stripping, d ynamical frictio n, and 
merging onto the large galaxy, can occur nea r luminous galaxies ( Kravtsov et al. 2004b: Bullock fc JohnstonI 



2005 



Zentner et aTll2005l ; [Maver et al.ll2006h . When a dwarf galaxy enters the environment of a large galaxy, 



tidal stripping can reduce its mass. In addition, the larger dwarf galaxies are preferentially dragged to the 
center and merge with the large galaxy. All of these processes alter the predicted mass function, but to 
follow them all requires difhcult-to-execute and physically uncertain simulations. From the point of view 
of observations, the ram pressure stripping removes any existing neutral gas disks, th e component of dwarf 



galaxies known to extend furthest out into the dark matter halo ( Stoehr et al. 20021 ). With only the stars 



it is difficult to probe any dark halo that might still surround a Local Group dwarf. 

On the other hand, isolated dwarf galaxies are possibly much simpler systems. Without a large galaxy 
nearby, the physical processes described in the previous paragraph cannot occur, simplifying the prediction of 
their mass functioii . In addition, as it happens, isolated dwarf galaxies essentially always have i ntact HI disks 
(jCeha et al.ll2006ar ). allowing us to probe their masses out to large radii with relative ease (jSwaters et al 
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20021 ). The disadvantage of isolated dwarf galaxies is that they are far away, making them difficult to find 



at the lowest luminosities. 

This paper represents a first step at examining the low mass end of the mass function for isolated 



galax ies. We use a sample of dwarf galaxies selected from the Sloan Digital Sky Survey (SDSS; lYork et al 



20001 ). From the SDSS, we evaluate their number densities, and we estimate their maximum circular velocity 
using single-beam HI profiles taken using the Arecibo Observatory and the Green Bank Telescope. We 
compare these observations to theoretical predictions for the number density of isolated halos of the same 
circular velocity. By comparing circular velocities, we avoid the considerable theoretical and observational 
uncertainties in determining the total mass of halos and galaxies. While this approach represents a beginning, 
we discuss in the conclusions how we can improve our observational estimates of the masses and number 
densities of these objects, and also find lower mass galaxies. With a much larger sample of isolated low mass 
galaxies, our technique can provide a stringent and unique constraint on the power spectrum at small scales. 

In Section [21 we describe the optical and radio observations our results are based on. In Section [3l we 
describe the TuUy-Fisher relationship for isolated galaxies. In Section 01 we describe our theoretical models. 
In Section [SI we compare the observations to the theory. In Section [S] we explore the robustness of our 
results to our definition of "isolation" in this context. In Section [3 we examine the relationships among the 
baryonic, stellar, and total mass estimates in our sample. In Section [51 we describe how these results might 
be improved upon in the future. Finally, in Section [9l we summarize our results. 

For determining luminosity distances and other derived parameters from observations, we have assumed 
cosmological parameters fio — 0.3, flA = 0.7, and Hq = 100 h km s~^ Mpc"^. Where necessary, we have used 
h = 0.7; otherwise, we have left the dependen ce on h explicit . AH m agnitudes in this paper are ivT-corrected 



to rest-frame bandpasses using the method of iBlanton et al.l (|2003l ) and kcorrect v4_l_4, unless otherwise 



specified. Because of the small range of look-back times in our sample (a maximum of around 700 Myr), we 
do not evolution-correct any of our magnitudes. 



Observations 



2.1. Sloan Digital Sky Survey 

To evaluat e the luminosi t y func tion of dwarf galaxies, we use a modified version of the SDSS spectro- 
scopic catalog. Blanton et al. ( 2005h describ e our sample, which is a subsample of the New York University 



Value-Added Galaxy Catalog (NYU-VAGC; 



Blanton et al 



20051) . We ha ve up dated that cata l og fro m SDSS 



Data Release 2 to Data Release 4 (DR4; Udelman-McCarthv et allbood ). The lBlanton et all |2005h catalog 
represents a significant improvement over naively selecting galaxies from the SDSS catalog, which is not 
optimized for nearby, low surface brightness galaxies. It is selected with an optical flux limit of ~ 17.8. 
For each galaxy, the catalog provides the SDSS redshift, emission line measurements, multi-band photome - 
try, structural measurements and environment estimates (for more ca talog details s ee 



Blanton et al. 



20051) . 



Distances are estimated based on a model of the local velocity field (jWillick et al.lll997f ). Distance errors 
have been folded into error estimates of all distance-dependent quantities such as absolute magnitude and 
HI mass. 

For our purposes, this catalog suffers from one maj or selection effect, due to the difficulty of detecting low 
surface brightness galaxies in the optical. As shown in Blanton et al. ( 20051 ). t he completeness as a function 
of half-light surface brightness drops below 50% at /i5o,r ~ 23.5 mag arcsec"^. [Blanton et al.l (|2005r ) present 
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a simple model for the effect that surface brightness has on the completeness, which assumes a log-normal 
surface brightness distribution with a mean that decreases as luminosity decreases. Figure [T] shows what this 
model implies about the missing fraction of galaxies as a function of luminosity in our sample, in terms of the 
correction factor c we must apply to recover the "correct" number densit y at each lum inosity. For comparison, 
we also show the equivalent factor for dwarf g alaxies in the Local Group ( Mateol l l998h giv en the SPSS surface 



bright ness cutoff, and in the local catalog of iKarachentsev et al.l (|2004 ). Note that the iKarachentsev et al 



(|2004[ ) curve is complicated by the fact that the quantities in the catalog (Holmberg radius and flux) are 
not enough to infer a half-light surface brightness (the terms in whic h we have calc ulate d the completeness) 
uniqu ely in general, even if we assume an exponential profile. The Mated (1998) and Karachentsev et al 
(|2004l ) catalogs bracket our correction factor o f 1.7 at Mr — S logm ~ 15, indicating it is roughly correct 
to within 20%. Therefore, we use the model of Blanton et al. ( 20051) to correct the luminosity functions we 
present here — without any a posteriori corrections. Nevertheless, the uncertainty of how many galaxies we 
are missing due to this effect is the most worrying one affecting our results. 



2.2. Environments of galaxies 

As argued above, there are advantages to finding isolated galaxies with which to test the circular velocity 
and mass functions. We cannot rely on the SDSS alone to determine whether a galaxy is isolated, for several 
reasons. First, the angular distances between nearest neighbor galaxies can be large for this nearby sample — 
for example, searching a 1 Mpc region around a galaxy 30 Mpc away corresponds to 2 degrees on the sky. 
Many of our dwarf galaxies are on the SDSS Southern stripes, which are only 2.5 degrees wide. In addition, 
because the SDSS reduction software is not optimized for large, extended objects and fails to process them 
correctly, the SDSS catalog docs not contain many of the bright galaxies within 30 Mpc. Thus, to calculate 
the environments of our dwarf galaxy sample, we need a supplemental catalog that extends beyond the SDSS 
area and contains the brightest galaxies. 

Both of these conside rations drive us to use the The Third Reference Catalog of Galaxies (RC3; 



de Vaucouleurs et al.Ml99ll ). which is a nearly complete catalog of nearby galaxies. To determine envi- 
ronments for our dwarf galaxies, we must determine the distance of each to its nearest "luminous" neighbor. 
In this context, we define galaxies as luminous when Mr — Slogj^g h < —19, corresponding to circular veloc- 
ities of Vc > 140 km s^^ for galaxies on the TuUy-Fisher relationship (see Section [3]). From the B and V 
photometry listed in RC3, we infer Mr for each galaxy. For galaxies which have the relevant entries listed, 
we call galaxies luminous if Mr — 51og]^o < —19. For galaxies which do not have the relevant entries, but 
do have HI data listed, we call them luminous if W20 > 300 km s~^ (as described in 5 12.31 W20 is twice the 
maximum circular velocity of the HI gas). Finally, there are some galaxies with neither HI data nor optical 
photometry listed in RC3. For this small set, we extract the "magnitude" from NED (which empirically 
is very similar to the B band RC3 magnitude for galaxies which have both) and guess Mr based on that 
magnitude. We call these galaxies luminous if Mr — 51og]^o ^ < —19. Additionally, we update the coordinates 
in RC3 using the NASA Extragalactic Database (NED) coordinates for each of the catalog objects. This set 
of bright galaxies is not perfectly uniform, but is suitable for our purposes. 

We combine the SDSS galaxies with Af,. — Slog^p/i < — 19 with the RC3 luminous galaxy catalog 
(removing repeats between the two). Then, we determine the nearest neighbor distance by asking whether 
there is a luminous RC3 or SDSS neighbor within 2 h^^ Mpc and 400km s~^ in redshift for each galaxy in 
our sample. In this paper, we will generally refer to isolated galaxies as those with no such neighbor within 
Him = 1-0 h^^ Mpc. However, in order to test the robustness of our results, we will vary our procedure 
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below by using alternate limits of Mr — 5 logj^Q ft, < —18 (a "fainter tracer sample" ) and Mr — 5 logj^Q h < —20 
(a "brighter tracer sample"). We will also vary the rp limit, as we describe explicitly below. 



We can calculate the number density of galaxies by weighting with the 1/Vmax values that lBlanton et al 



(|2005l ) describes. Of course, some of the volume actually is not available (e.g. in a cluster a galaxy will 
never appear isolated) and in principle we might need to account for this fact. To do so, we can measure 
the "isolated fraction," the fraction of the volume of our sample which is isolated from luminous galaxies, by 
randomly placing points within the volume and testing their environment. This fraction is 0.95, and varies 
by less than 0.05 when we vary our definition of "isolated" as noted above, minor compared to our other 
uncertainties. Thus, we simply ignore this correction, and calculate the number density of isolated galaxies 
over the entire volume. 

The thick line Figure [2] shows the cumulative number density of such isolated galaxies as a function of 
absolute magnitude, corrected for surface brightness incompleteness. This statistic depends on the definition 
of "isolated" that we use; the two thin lines use the fainter and brighter tracer samples described above. As 
this comparison shows, a change in magnitude of the tracer by about 1 mag is equivalent to a 20% change 
in the cumulative number density at Mr — 5 log^Q h ^ —15. 

Taken together these data show that the number density of isolated galaxies with — 5 log^^Q h < —14.7 
is (6.83 ± 2.3) X 10~^ Mpc~'^, with the uncertainty dominated by the definition of "isolated" and the 
fraction of galaxies missing due to surface brightness effects. 



2.3. Radio observations at 21cm 



Geha et al.l (|2006a[ ) present the radio observations of dwarf galaxies selected from the SDSS sample 
described above. The data were obtained on the Green Bank 100-m Telescope (GET) and at Arecibo 
Observatory. Each observation had a velocity resolution better than about 3 km s~^. The optical half-light 
radii of the dwarf galaxies in our sample are typically ~ 8", and should be completely contained with the 
radio beamsize of 3' and 9' for Arecibo and GET, respectively. Here we will restrict ourselves to galaxies with 
detected HI (which iGeha et al.ll2006al showed dominated the isolated galaxy population) and with b/a < 0.5 
to minimize inclination effects. 

We compute the 20% HI line- width (W20) by finding the peak HI flux within 150 km s~^ of the optical 
radial velocity of each galaxy and computing the difference between the nearest points having 20% of the 
peak flux. The integrated HI flux is calculated by expanding the W20 values by 20 km s^^ on each side and 
integrating the flux in this region. Errors bars on the line widths and integrated fluxes were computed using 
a Monte Carlo bootstrap method: noise was added to the stacked one-dimensional radio spectra (based on 
the observed variance in the baseline) and the observed quantities remeasured. We calculated error bars 
on the line-width and integrated flux from the scatter in the mean quantities recovered from Monte Carlo 
simulations. 



As lCeha et al.l (j2006a[ ) showed, the dwarf galaxies in our sample have a signiflcant rotation component. 
We derive the maximum rotation speeds as follows. We correct the observed HI line- widths for line broadening 



due to turbulent velocity dispersion and inclination using the formula first proposed by iBottinelli et al 
(Il983h : 

W20 - W20t 



(1) 



where W20 is the observed HI line-width, W20t is the turbulent velocity correction term and i is the 
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inclination angle inferred from the optical images. We confirm the validity of a line ar turbulence cor rection 
by modeling the integrated velocity profiles of simulated galaxies constructed from 'HernquistI (Il993[) model 
disk galaxies. For nearby dwarf galaxies with rotation velocities similar to our sample. Begum et al.l 12006) 
have measured a velocity dispersion in the gas component of aios = 8 km s^^ from 2D velocity maps. Using 
the Begum et al. value, this results in a turbulence correction of W20t = 16 km s~^, which we use here. 
Altering this turbulence correction (say to 25 km s~^) does not change our results below. 

Figure [3] shows our best estimate of Knax for our dwarf galaxies (with b/a < 0.5) as a function of the 
distance to the nearest luminous luminous galaxy. Clearly there is a strong relationship between VJnax and 
projected separation. Within 1 h^^ Mpc of luminous galaxies there is a population of objects with Vmax < 40 
km s^^, that is much rarer at large separations. Our galaxies are selected to be distributed roughly evenly 
in the range —13.5 > Mr — 5\ogiQ h > —15.5, independent of environment. Thus, at these luminosities the 
"forward" TuUy-Fisher relationship — the circular velocity at a fixed luminosity — appears to be a strong 
function of nearest neighbor distance. This result does not imply that no low circular velocity galaxies exist 
in the field, simply that any such galaxies are too low luminosity to make it into our sample. 



As iGeha et al.l (j2006aj) and others have found, dwarf galaxies near a luminous neighbor also tend to be 
red and gas poor. Taken all together, these results suggest that some important physical effects are shaping 
the gas content, star-formation histories and inferred dynamics of dwarf satellite galaxies relative to isolated 
dwarfs. For this reason, we choose to concentrate our attention here on the isolated galaxies, whose Vmax 
values have a smaller dispersion, and whose properties in general we expect to be less altered since formation, 
relative to those dwarfs perturbed by a massive neighbor. 

Although it is immaterial to our analysis below, it is interesting to ask what physical effects are causing 
the trend in Figure [3l We can think of three explanations. First, for dwarfs in the vicinity of luminous 
galaxies, ram pressure stripping could remove gas at the largest galactocentric radii first, reducing the 
maximum circular velocity traced by HI emission. Second, tidal stripping could reduce the total mass and 
thus the circular velocities. Third, interaction-triggered star- format ion could either raise the luminosities 
of satellite galaxies relative to isolated galaxies, bringing lower mass systems into our sample if they are 
near bright galaxies, or speed up star- formation and use up the gas in the outer disk of the galaxy. To 
help investigate this question, in Figure [3] we have distinguished between sin gle-peaked HI profi les (open 



circles) and double-peaked or flat-topped profiles (filled circles) , as classified by iGeha et al.ll2006a[ ) . Because 
the dwarfs with low Vi„ax are predominantly single-peaked profiles (and we have enough resolution to see 
double-peaked profiles if they existed), we favor the explanation that gas has been stripped from the outsides. 
However, a definitive conclusion awaits a comprehensive analysis, including more detailed dynamics of these 
galaxies. 



Tully-Fisher relation for isolated galaxies 



Here we give an estimate of the Tully-Fisher relationship for isolated galaxies. We base our estimate 
on the sample descri bed i n the previous section for low luminosities, plus the recently published samples of 
Pizagno et al. ( 2006h and Springob et al. ( 2005 ) for higher luminosities. 



Pizagno et al.l (j2006f ) presented a Tully-Fisher survey of galaxies found in the SDSS, using follow-up Ha 
rotation curves. They have performed model fits to the disk components of these galaxies to inch nation- 
correct their maximum circular velocities. Unlike previous Tully-Fisher samples (e.g., Courteaulll997 ). which 
were selected to be very homogeneous sets of galaxies, the Pizagno et al. ( 2006h sample spans a large range 
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of galaxy types and colors, resulting in a somewhat larger scatter in the Tully-Fisher relationship than found 
in other studies. For each of their galaxies we have determined its environment in th e same manner as for 
our low luminosity sample, and we only consider isolated galaxies here. However, as iPizagno et al. 
show, and we have confirmed using our own measurements of environment, there is very little dependence 
of the Tully-Fisher relation on environment at high luminosities. As with the low lumi nosity sample, we 
restrict ourselves to galaxies with b/a < 0.5, leaving 35 galaxes from Pizagno et al. ( 2006f ). 



Springob et al.l (|2005f ) have compiled HI spectra from archival data sets for around 9000 galaxies in 
the local Universe, observed originally with Arecibo, the 91m and 42m Green Bank telescopes, the Nangay 
telescope, and the Effelsberg 100m telescope. They have homogeneously analyzed these spectra, measuring 
their widths and HI fluxes. All their galaxies have optical data associated with them, including a measure 
of their axis ratio (b/a). We use the turbulence and inclination corrections of Equation ([1]) to convert the 
W20 values to Vmax. In addition, in order to get luminosities and environments for these galaxies, we match 
this list to the full low-redshift catalog from SDSS DR4 (just as for our dw arf sample) . Finally, w e restrict 
this sample to the isolated galaxies with b/a < 0.5, leaving 85 galaxies from ISpringob et al. I (|2005h . 



Figure m shows the Tully-Fisher relationship for all three samples of galaxies. The low luminosity points 
are from our isolated sample. The crosses in Figure d] show the median Vi„ax values in several bins centered 
on Mr — blogiQh = —14.7,-18.5, —19.5, and —20.5, and 1 magnitude in width. Table [T] lists the median 
velocity in each bin (and several other bins needed below). 



of 



4. Theory 

In order to find an ap propriate theoretical comparison, we use the iV-body, pure dark matter simulations 



Kravtsov et al.l (|2004al) . They simulated a cubic box 80 Mpc on a side, with 512'^ particles, each 
around 3.2 x 10* Mq in mass. The world model and transfer function for the simulation correspond to 
flm = 0.24, = 0.76, rib — 0, h = 0.73, n = 0.95, and as — 0.75. Although this cosmology is not precisely 
the current best-fit model, we describe below how we implement small corrections to account for this fact. 
We used the outputs from z = 0. 



Dark matter halos were identified using the method described in iKravtsov et al.l (j2004al ) . They calcu- 



lated the maximum circular velocity of each halo by determining the enclosed mass M{< R) as a function 
of radius, using the spherically symmetric approximation v = ^ GM (< R)/ R, and determining the peak in 
the rotation curve. Strictly speaking, this T^nax value is only comparable to our observed galaxies if their 
21cm emission probes this peak. However, because most of our galaxies show flat-topped or double-peaked 
HI profiles, we are probably close to satisfying this condition. 

To define "isolated" we use tracer halos with Vniax values corresponding to the Tully-Fisher results 
listed in Table [T] We project the distribution of halos in a random direction, and perturb the "redshift- 
space" positions of the halos to account for their peculiar velocities. Then we define "isolated" in the 
simulation using the same geometrical considerations used in the observations (|Ai;| < 400 km s~^ and 
Tp < 1 h^^ Mpc) but relative to the tracer halo population. The velocity function of these isolated halos 
is shown in Figure [5] as the histograms. Down to 70 km s^^ each histogram comes from the simulation, 
but below that we extrapolate the velocity function as a power law (roughly <!>(> T^nax) oc V^^^). The 
three histograms, as labeled, correspond to tracers with VJnax > 110 km s~^ (corresponding to galaxies with 
Mr — 51og]^Q/i < —18), Knax > 140 km s^^ (corresponding to galaxies with Mr — blogi^h < —19), and 
VJnax > 180 km s^^ (corresponding to galaxies with Mr — 51og]^Q/i < —20). For the bulk of this paper. 



- 8 - 



we will be concerned with the central class, shown as the thick histogram, but will use the other results to 
quantify how much our results depend on the choice of tracer population. 

Because of the complex geometrical definition of "isolated," it is useful to have this iV-body estimate 
of the mass function. However, the cosmology used for the simulation does not correspond precisely to 



the current best fit cosmology (e.g. IXegmark et al.l 120061 ): in particular it does not include the effects of 



baryons on the initial power s pectrum. We adj ust for this difference by using the excursion set formalism 
and the transfer fu nctions of Eiscn stein &: Hul (jl998l ). along with the mass function approximation from 



Warren et al.l (|2006r ). These methods are able to predict the mass function for any hierarchical cosmology. 



as a function of redshift and of large-scale environment. We use a particular implementation provided by 
Andreas Berlind (private communication). These methods yield the mass function of galaxies, which we 



(2001 



3), using M* = 1.5 X lO^^ 



convert into a circular velocity function using the methods of iBuUock et al 
/j-i Mq and nm = 0.24. 

First, we need to evaluate what large-scale environment our definition of "isolated" corresponds to. 
The smooth line in Figure [5] is the prediction for th e mass function of halos in large-scale underdensities of 
5 = —0.4, for the cosmology used in lKravtsov et al.l (|2004a[) . From this agreement, we conclude that in the 
excursion set mass functions, S = —0.4 is the underdensity that is most comparable to our isolation criterion. 

Second, in order to adjust the results of^ Kravtsov et al. ( 2004a) for co smology, we evaluate the ratio fc 
between the ve l ocity function for the cosmology of iKravtsov et al.l (|2004a|) (listed above) and the best fit of 
Tegmark et al.l (|2006[ ) (the only differences are that in the latter erg = 0.76 and fibft,^ = 0.022). Figure [6] 



shows this ratio as a function of maximum circular velocity. Over the ra nge we will use here, th is correction 
is never more than about 10%. In order to compare our results from the IKravtsov et al.l (j2004al ) simulations 
to observations, we first apply this correction factor to the predictions from the simulations. 



In addition, we have incorporated a "warm dark matter" version of the Tegmark et al. ( 20061 ) cosmology 
that is identical in its large-scale structure, but includes a light dark matter particle, with uidm = 0.5 keV. 
The lightness of this particle increases its free streaming length, which smoot hs fluctuations on small scales. 
Here we apply the adjustment required to the transfer function as outlined by Abazaiian ( 20061 ). We simply 
input this new transfer function into the excursion set calculation of the mass function. Of course, on the 
scales compa rable to the free-s treaming length, the collapse of structure will cease to be hierarchical, as 
described bv iBode et al.l (|200lh . Thus, almost by deflnition this prediction will be incorrect; however, we 
use it as a rough approximation to a more correct prediction which might be available in the future. Figure 
|6| shows the "correction" as a function of Vmax as the dashed line. In order to predict the warm dark matter 
case, we apply this correction factor to the simulations. 

Figure |7| shows the resulting Vmax functions in the cold dark matter and warm dark matter cases. In 
the next section, we describe how we compare these theoretical predictions to the observations. 



5. Comparing theory to observations 

The simplest possible comparison we can make between simulations and the observations is to try to put 
some observed points on the velocity function. After all, we have measured the number density of galaxies 
as a function of luminosity, and from the TuUy-Fisher measurements, we know the relationship between 
luminosity and the circular velocities. Given the luminosity function of Figure |2| (using the tracers with 
Mr — 51ogiQ h < —19), plus the relationship between luminosity and Vmax from Figure HI we also plot the 
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observed number density as a function of Vmax as the four points. In this comparison, there appear to be 
shght discrepancies between the cold dark matter model and the observations, particular at the bright end. 
The warm dark matter model is nearly as good a fit to the data but somewhat underpredicts (at a bit more 
than 1(t) the number of isolated low circular velocity galaxies. 

However, this method of comparison is sensitive to bias related to scatter in the relationship between 
luminosity and circular velocity. A more robust comparison can be achieved as follows. For a given predicted 
circular velocity function and a given observed luminosity function, we can find the relationship between cir- 
cular velocity and luminosity that makes them consistent with one another. We do so here by parameterizing 
it as a piecewise linear relationship Afr(Vinax) between the circular velocity and a mean absolute magnitude, 
with Gaussian scatter about that mean. We vary the parameters of the piecewise linear relationship to fit 
the luminosity function (using the Levenberg-Marquardt method implemented in the IDL routine mpfit 
distributed by Craig B. Markwardt). However, we fix the Gaussian scatter to have an = 2.5 for Vmax < 10 
km s~^, to have ctm = 0.4 for Vmax > 100 km s^^, and to vary linearly with circular velocity in between. 

Figure [5] shows the conditional distribution of V^max a-s a function of M^, given our best fit relationship 
(using tracers in the simulation with Vljiax > 155 km s^^ and corresponding tracers in the observations 
with Mr — 5\ogiQ h < —19). The lines are the quartiles of the distribution. The overplotted points are the 
data from Figured Clearly the median values (shown as the boxes) agree rather well with the predictions. 
This constitutes a confirmation that the TuUy-Fisher relationship and the luminosity function together are 
consistent with cold dark matter predictions. 

How well can these data exclude alternate scenarios? We explore this question by considering the warm 
dark matter model described above, with a dark matter particle mass of 0.5 keV. We perform the same 
procedure as described above and obtain the predictions shown in Figure [9] Here, the observed circular 
velocities of low mass galaxies tend to be higher than predicted, but not by significant amounts. The median 
is about 4(7 away from the predicted median, which given the systematic uncertainties here we regard as a 
marginal exclusion of this model. 

6. Robustness relative to our definition of "isolated" 

Here we examine the robustness of our results to our definition of "isolation," and how we relate isolation 
in the observations to isolation in the simulations. There are a number of arbitrary decisions we have made 
here to define "isolated" in the observations. In particular, we chose a certain projected distance rp from 
galaxies of a certain absolute magnitude Af^.bright (or brighter). We must examine the sensitivity of our 
results to these arbitrary choices. In addition, we have also had to define "isolated" in the theoretical 
predictions. To do so, we have had to relate the absolute magnitude Mr used above to a T4iax of the halos 
in the simulation. Of course, we do not know the exact correspondence, and so we need to examine how our 
results depend on errors in our estimate of it. 

Figure [TOl examines the sensitivity of our results to both sets of decisions. The top panels show the ratio 
of the predicted <f>iso(> Knax) at Vi„ax = 56 km s~^, to the observed <I>iso(< Mr) at Mr — 51ogio h — —14.7. 
The bottom panels show another comparison of the observations to the theory: the ratio of the observed 
Vmax = 56 km at Mr — 51og]^Q h = —14.7 to that predicted by equating the number density of halos 
larger than a given VJnax to the number density of galaxies brighter than Mr — 51og]^Q h — —14.7. The left 
panels correspond to our cold dark matter model, and the right panels correspond to our warm dark matter 
model with to_dm = 0.5 keV. 
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There are 27 points shown on each plot, each corresponding to a different definition of "isolated" in the 
observations and the theory. First, we check three different choices of tracer sample (Afr,bright — Slogj^g ^ — 
— 18, —19, and —20), as shown by the rough horizontal position. Second, for each of these three choices of 
galaxy sample, we choose three different Knax values with which to define our halo sample for comparison. 
Our choices of predicted circular velocities for each observed sample are listed in Table [TJ In Figure [TOl the 
larger symbols correspond to higher circular velocities. Third, for each of those nine choices, we choose three 
different choices of rp (0.7, 1.0, and 1.3 Mpc), using the same radius for theory and observation. These 
three cases are offset from each other slightly in Figure [10] (left to right, respectively) for clarity. 

From these results, it is clear that our systematic uncertainties are (fractionally) about 0.3 in the 
comparison of number densities, and 0.1 in the comparison of circular velocities. The difference is at least 
partly due to the approximate dependence $(> Vmax) cx V~^^'^ . In addition, we expect the comparison of 
circular velocities to be more robust, since it depends less on the scatter in the relationship between luminosity 
and circular velocity. While for our sample, these systematics are about equivalent to the systematics 
associated with our surface brightness completeness selection, the systematics shown in Figure [TOl illustrate 
what will ultimately be the most difficult uncertainty to overcome, even when much more complete galaxy 
samples are available. 



Baryonic and stellar content of low mass galaxies 



In the previous section, we showed that the cold dark matter model reasonably explains the number 
densities and circular velocities of low luminosity galaxies. Assuming that the relationship between circular 
velocity and total mass that cold dark matter theory predicts is correct, we can now investigate the baryonic 
and stellar mass content (relative to the total mass) of these low mass galaxies. This census of the matter 
in dwarf galaxies may help us understand their creation and development over time. 



To infer the total mass from Knax we use the methods of iBullock et al.l (j2001bl ). using = 1.5 x lO'^^ 
h^^ Mq and fim = 0.24. These methods have been calibrated down to masses of 10^"'^ Mq, or about 
85 km s~^, so our use of them for the lowest luminosity galaxies represents an extrapolation of the current 
theoretical understanding. At our typical Knax ^ 56 km s~^, the virial mass of the halo determined by this 
method is 2.5 x 10^" h-^ M^. 



T o infer the stellar mass, we use the optical broadband SDSS data and the methods of Blanton &: Rowei^ 
([20071). Any method for inferring the total stellar mass is sensitive to the initial mass function (IMF) of 
stars, since the lowest mass stars (most of those below 0.5 Mq) contribute almost no optical light but are 
a sig nificant fraction of the ma ss. The differenc es between differ ent reasonable choices of IMF can be up to 
50%. iBlanton fc RoweisI ([20071) have chosen thelChabrieii (120031) IMF , and find stellar masses within about 
30% of those found for the same galaxies bv iKauffmann et al.l (j2003( ). using spectroscopic techniques. The 
median stellar mass for our isolated dwarf galaxy sample is 2.2 X 10^ h-^ Mq. 

By "baryonic mass," we mean here the sum of the stellar mass and neutral gas content, which we take 
to be Mf, — A/» + 1.4Mhi, where Mhi is the neutral hydrogen mass inferred from 21cm observations and 
the factor 1.4 accounts for helium, molecular clouds and metals. The median baryonic mass so defined is 
2.5 X 10^ Mq — much larger than the stellar mass contribution. Naturally, there may also be ionized 
hydrogen in the galaxy, which we do not try to account for here. Note for the data from lGeha et al. I liOOGa) 
we do not try to correct for self- absorption of HI, because l ittle evidence for any in clination dependence of the 
HI to stellar mass ratio is found in our sample. However, ISpringob et al.l (|2005[ ) did make such corrections. 
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which can be up to 20%. 



For the samples of isolated galaxies used in this paper, Figure [TT] shows the ratio of baryonic mass to 
total mass as a functi on of r-band absolute magnitude (for h = 0.7). The dashed line is the cosmic mean 
b ased on the resul t s of Tegmark et al.l (j2006l ) (rifj/fim = 0.17) and the dotted line is t he mean for the sam ple 
of ISpringob et al.l (|2005f) taken alone. The mean of the low luminosity galaxies from lGeha et al.l (|2006al) is 
somewhat less than that of higher luminosity galaxies, about 8% of the cosmic value rather than 14%. The 
difference in the treatment of self-absorption may account for about half of this difference. 

The baryonic fraction may continue to decrease at lower masses. However, at Knax ~ 50 km s"-'^ iso- 
lated low luminosity galaxies do not show much evidence that they have expelled or ionized very much more 
of their cold gas than have their more massive counterparts. This measurement supports detailed models 
of the physics of gas blow-out and blow- way due to supernoyae, which predict a loss of only a few per - 
cent for galaxies in this mass range |Mac Low fc Ferraralll999t iFerrara fc Tolstovlbood Istinson et aPboOTi) . 



onic mass loss (Dekel & Silk 1986; Cole et al. 2000; Mori & Burkert 


2000; 


Bullock et al. 


2000; 


Benson et al. 


2003; 


Dekel & Woo 


2003; Tremonti et al. 2004; Croton et al. 2006). 


Of course, reasonable modifications to 



those models in which feedback prevents the gas from forming stars but keeps it mostly in neutral form 
and within the galaxy disk are probably tenable. Fu rthermore, althou gh some investigators have invoked 
outflows to explain the mass-metallicity relationship, iDalcantonI (j2007l ) have shown that alternate models 
without significant outflow can explain the observations. 

We can also look at the relationship between the stellar mass of the galaxies and their total mass. 
Figure [T^] shows the dynamical to stellar mass ratio as a function of stellar mass for the three samples used 
here. This relationship shows a strong trend, illustrating the strong dependence of star-formation efhciency 
on mass, at least for disk galaxies. It is, of course, exactly this dependence that causes the luminosity and 
stellar mass functions to be shallow while the total mass function of galaxies is so steep, as described in 
Section [TJ 



8. Future directions 

The analysis of this paper, while consistent with the CDM model, puts only mild constraints on alterna- 
tive models (for example, a warm dark matter particle as light as 0.5 keV is barely ruled out). How can this 
analysis be improved in the future? Two paths are possible: first, increasing the depth and completeness of 
our optical plus HI sample; second, using upcoming blind HI surveys to push to considerably lower masses. 
Both paths require improving our understanding of the theoretical predictions at the low mass end. 

Our analysis of the current SPSS s ample would be substantially improved with more HI follow-up 



observations. Because lGeha et al.l (|2006aj ) were studying the general HI properties of dwarfs, we only targeted 
about 12 systems edge-on enough and isolated enough to include in the analysis of this paper. Our results 
could be put on a much firmer footing with an increase in our follow up sample. In DR4 there are 64 
galaxies with b/a < 0.5, Mr — 51og]^o ^ > ~15 and rp > 1 Mpc, and obviously there are still more in 
later releases. Additional follow-up or deeper HI sky surveys should fill this gap in the future. 

The optical SDSS sample on which our analysis here is based will increase by DR8, the final SDSS 
release, perhaps by a factor of two. However, this will not substantially reduce the uncertainties in the 
luminosity function, which are already dominated by the surface brightness completeness correction (Figure 
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[T]). iBlanton et all |2005h concluded, based on introducing simulated galaxies into raw SDSS imaging data, 
that the SDSS photometric pipeline (optimized for reasonably high surface brightness galaxies around z ~ 
0.1) was probably failing at a brighter surface brightness limit than the data required, and that with a 
differently optimized pipeline it might be possible to push the surface brightness limits to 25 mag or more. 
Doing so would allow us to probe magnitudes as faint as My. — Slogj^g ^ ^ ~12 over cosmological volumes. 
Another p ossibility would be to wait for upcoming, deeper surveys s uch as the Dark Ene rgy Survey (DES; 
Wester et al. 2005), which has first hght in late 2010, Pan-STARRS 4 |Hodapp et al.l|2004t ), whose prototype 



telescope PSl will probably see first light in 2007, or the Large Scale Synoptic Telescope, which has first 
light in 2014. Searches for low surface brightness galaxies (through their diffuse light) tend to be dominated 
by the scattered light background, so it is difficult to anticipate how well of any of these surveys can do. 
Any of these possibilities would require spectroscopic follow-up — probably searching for 21cm emission in 
the radio or targeting HII regions in the galaxies. 

Blind searches for galaxies in 21 cm may show even more promise, particularly since isolated dwarf 
galaxies virtually always exhibit HI (jGeha et al.l l2006al ) and since the dynamics of each galaxy will be 
measured simultaneously with its detection. Unfortunately, HIPASS appears to not be deep enough to 
provide a competitive samp le in this resp ect (many of the SDSS galaxies in our sample are undetectable in 
HIPASS; ICeha et al.ll2006al ). ALFALFA iGiovaneUi et al.ll2005[ l can detect galaxies with HI masses of lO"^ 

M© at distances of 20 hr^ Mpc. If the full 7000 deg^ planned survey is completed, the overall volume 
mapped will be about 6 x 10'^ /i'^ Mpc"'^. That large a volume has an expected cosmic variance of around 
30%, though given our restriction to isolated regions the actual cosmic variance uncertainties will be lower 



Assumi ng a baryonic to total mass fraction of 0.02 (see Figure [TT]) and using the methods of iBullock et al 
(j2001a[ ). t his mass corresponds to Vmax ~ 20 km s^^. Galaxies of this mass in the preliminary ALFALFA 
catalog of iGiovanelli et al.l (|2007l) have a median W50 measurement of ~ 40 km s~^, consistent with this 



estimate. Assuming (conservatively) that 30% cosmic variance errors dominate the uncertainties, using our 
techniques here one could marginally exclude a warm dark matter particle of 2 keV in mass. Future surveys 
such as MIRANDA will contribute a similar volume (of a distinct chunk of the Universe), and increase this 
precision somewhat. In any case, these new HI surveys will push this technique into a low circular velocity 
regime that is currently only tested with observations of Local Group satellites. 

Even to make use of the current data, however, we probably require a better understanding of the 
theoretical predictions. For example, as we noted above, the excursion set predictions we are making for 
"warm dark matter" are not entirely self-consistent, since we expect in this regime that the hierarchical 
picture will start to break down. Under these conditions, iBode et al.l (j200l|) found that the number of 
forming halos was much smaller than the excursion set prediction. Thus, it may be that our observations 
would put stronger constraints on a correctly calculated prediction. However, doing so is difficult, since the 
effective mass resoluti on for warm dark ma tter simulations appears to scale much less favorably than for 
cold dark matter (see Wang fc Whitel 2007 , who indeed argue that even the number of halos predicted by 
Bode et al.ll200ll is an overestimate). 



Finally, it is worth noting two possible fundamental limits to the technique we describe here. First, we 
rely on the 21cm emission to probe the dynamics in the flat part of the rotation curve. While this appears to 
be the case for most of our galaxies here (based on their double-peaked morphology) , it will not necessarily 
be true of the typical very low mass galaxy. If so, one would need to resort to comparing to cold dark 
matter predictions in the very inner parts of halos. Second, at lower masses, it may happen that reionization 
evaporates the lowest mass halos, preventing the formation of stars and the existence of any neutral gas at 
all. If so, our counting technique will fail. For the galaxies in the mass range we study here, this appears to 
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be a weak effect, but it may occur at smaller masses. 



9. Summary 

We demonstrate that the predictions of cold dark matter are consistent with the number density of 
isolated low circular velocity galaxies (Knax ^ 50 km s~^), at a precision of about 30% (Figures [7] and 
[8]). Our major systematic uncertainties (which dominate our error budget) are related to our definition of 
"isolated" and our model of the surface brightness completeness of the SDSS at low luminosities. These 
results represent a valuable independent check of the small-scale predictions of the cold dark matter model. 
Our technique avoids many of the systematic uncertainties associated with observations of satellite galaxies 
in the Local Group, related to the complex physics of ram pressure and tidal stripping. From a statistical 
point of view, our results are less powerfu l than th e current constraints from the Lya forest power spectrum 



([Narayanan et al.l bOOOt lAbazaiianl 120061 : ISeliak et al. 2006,) , and only marginally exclude a dark matter 



particle with moM ~ 0.5 keV. With an improved sample, this technique will offer a powerful and unique 
constraint at small scales. 

We also find several secondary results that are relevant to the formation of these dwarf galaxies: 

1. At low luminosities, the TuUy-Fisher relationship appears to be a function of environment, with dwarf 
satellite galaxies having lower circular velocities than isolated dwarf galaxies, by up to a factor of two 
or more. There is circumstantial evidence from the nature of the HI profiles that this effect is due to 
stripping of the outer gas in satellite galaxies, though other processes may be at work. 

2. The baryonic mass fraction of galaxies (counting neutral gas plus stars) appears to be a weak function of 
luminosity down to Mr — 51og]^o ^ ~ ~14, decreasing by 40% at most (from about 14% to at minimum 
about 8% the cosmic mean; Figure [T2| . This result disfavors models which call for a preferentially 
large amount of baryonic outflow in dwarf galaxies (due to internal processes such as feedback). 

3. The ratio of total to stellar mass is a very strong function of stellar mass, ranging from 50 or so at the 
highest luminosities to over 1000 at the lowest luminosities. 



Taken together with the deficit of neutral gas in dwarf galaxies near luminous neighbors (jGeha et al.ll2006a^ . 
these results suggest that the primary effects that remove gas and end star-formation in dwarf galaxies (with 
circular velocites of 50 km s"""^ or so) are external interactions with bright neighbors, rather than internal 
processes such feedback and outflows. 

Increasing our sample of isolated dwarf galaxies with HI follow-up (it is straightforward to increase the 
current sample by more than a factor of five) would improve the precision of all of these results, which are 
based on a relatively small sample of isolated galaxies. Upcoming blind HI surveys such as the ALFALFA 
survey on Arecibo are going to impose better constraints on cosmological models at small scales, as well as 
better explore the issues of dwarf galaxy formation. They will propel the study of field dwarf galaxies into 
a low circular velocity and low mass regime previously studied only in the Local Group. 
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Fig. 1. — Estimates of necessary correction factors c f or surface brightness incompleteness as a function of 
absolute magnitude. The solid line is the estimate from Blanton et al. ( 2005h using a simple model of surface 
brightness as a function of absolute magnitude and determinations of the completeness of SDSS as a function 
of surface brightness. The dashed li ne uses the distribution of surface brightnesses as a function of absolute 
magnitude from the local cat alog of Karachentsev et al.l (j2004l ). The dotted line uses the distribution from 



the Local Group according to iMateol (|l998r ). 
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Fig. 2. — Cumulative luminosity function of isolated galaxies, as defined in the text. The meaning of 
"isolated" depends on the tracer galaxies used. The thick histogram represents the median relationship 
using tracer galaxies with M,. — Slogj^Q h < —19. The thin histogram represents a change of 1 magnitude in 
the absolute magnitude limit used for the tracer galaxies; the upper histogram uses brighter galaxies, and 
the lower histogram uses fainter galaxies. 
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Fig. 3. — Estimate of VJnax for dwarf galaxies in our sample as a function of distance to the nearest luminous 
neighbor galaxy. Filled symbols are galaxies with flat or double-peaked HI profiles, open symbols are galaxies 
with single-peaked HI profiles. 



- 21 - 



o 



2.8 
2.6 
2.4 
2.2 



Of) 

o 2.0 



1.8 
1.6 



— I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — p 
• Geha et al. (2006) 

o Pizagno et al. (2006) 
Springob et at. (2005) 

+ Binned median values (boundaries shown) 



r o 



1.4 



_l I I I I I I I I I I I I I I I I I L_ 



12 



14 



16 



18 



20 



22 



Fi g. 4. — TuUy-F isher relation for isolated galaxies. Large open symbols a re from the optical rota tion curves 

(|2005h . and filled 
max values in several 



of iPizagno et al. I 12006 ). small open s ymbols are from th e compilation of ISpringob et al 



symbols are from the HI linewidths of iGeha et al.l (j2006bf ). Crosses show the median Vr 
bins; the bin centers and median values are listed in Table [T] 
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Fig. 5. — Cumulative Knax function of isolated halos from A^-body simulations, for various choices of tracers. 
The thick histogram is the abundance of halos in the simulation subject to the isolation criterion described 
in the text, based on the projected distance and redshift different relative to the nearest, other "tracer" 
halo, for tracers with Vmax > 140 km s^^. The thin histograms explore the effect of trying d ifferent tracer 



populations, as labeled. The smooth line is an approximation based on the mass functions of I Warren et al 
(|2006f ). described more fully in the text, using a relative large-scale bias oi 5 = —0.4 for the halos. 
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Fig. 6. — Correction to the mass function between the cosmology of the A^-body simulations of lKravtsov et al 



(|2004al ) and the cosmologica l mod els we want to test. The solid line is the correction to the standard 
cosmology of Tegmark et al. ( 2006[l based on WMAP and SDSS large-scale data. The dashed line is the 
correction to that same cosmology, but with a light dark matter particle (todm = 0.5 keV). Correctio ns are 
based on the ratio between the appropriate cosmologies in the approximations of IWarren et al.l (|2006l ). 
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Fig. 7. — Cumulative Vmax function of isolated halos from Figure [5l corrected to the iTegmark et alj ( 20061 ) 
cosmology (upper histogram) and that same cosmology with a light dark matter particle (lower histogram) . 
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Fig. 8. — Conditional distribution of circular velocity as a function of absolute magnitude, to make our 
observed isolated galaxy luminosity function consistent with the predicted isolated halo circular velocity 
function. Lines are the quartiles of the distribution. The overplotted points and binned median values are 
the data from Figure ID 
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Fig. 9. — Similar to Figure [8] but based on the warm dark matter models. 
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Fig. 10. — Dependence of our results on our definition of "isolated." The left panels refer to our cold dark 
matter model. The right panels refer to the warm dark matter model with moM = 0.5 keV. The top 
panels show the ratio of the observed number density of isolated galaxies with Mr — Slogj^p ^ < —14.7 to 
the predicted number of isolated halos with Vmax > 56 km s^^ (that is, the vertical offset in Figure[7]). The 
bottom panels show the ratio between the median observed circular velocity at M^- — hlog^Q h < —14.7 to 
that predicted in Figures [8] and [Q] All of these results are shown as a function of how we define "isolated," 
which we do in 27 different ways here. We choose: (a) three different absolute magnitude limits for the 
tracers in the observed sample, Mr — Blog^gh < —18, —19, and —20, as shown by the rough horizontal 
position; (b) for each sort of observed tracer, three choices of minimum Knax for tracers in the predicted 
sample, as listed in Table [1] and as shown by the size of the symbols (larger corresponds to higher circular 
velocity of tracer); and (c) three different minimum projected radii from the tracers, rp = 0.7, 1.0, and 1.3 
Mpc, as shown by small offsets in the horizontal positions (left to right). 
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Fig. 11. — "Baryonic" mass relative to total mass as a function of luminosity for the samples of isolated 
galaxies with b/a < 0.5 in the two HI samples used here, as marked. Baryonic mass is defined as the 
stellar mass plus the neutral gas mass, as described in the t ext. The dashed line at 0.17 is the cosmic mean 
based on cosm o logica l measurements ( Tegmark et al. 20061) . The dotted line at 0.025 is the mean of the 
Springob et al. ( 2005h measurements. 
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Fig. 12. — Dynamical to stellar mass ratio as a function of stellar mass, for the three samples used in this 
paper (restricting to isolated galaxies with b/a < 0.5). 
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Table 1. Comparable circular velocities and absolute magnitudes 



Mr — 5 logio h 


V^nax (km S ^) 


Vmax (km s ^) 




(from TF) 


(used) 


-14.7 


56 ±3 




-18.0 


108 ±5 


95, 110, 125 


-18.5 


116 ±6 




-19.0 


143 ±7 


125, 140, 155 


-19.5 


156 ±7 




-20.0 


180 ±4 


170, 180, 195 


-20.5 


214 ±7 





Note. — For each choice of Mr (first column), 
this table yields the comparable Knax resulting from 
matching to the TuUy-Fisher relation in Figure [3] 
(second column, "TF" ) . It also lists the V^nax values 
we actually tried in Figure [TO] (third column). 



